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Electronic transport through a 2D deca-nanometer length channel with correlated and anti- 
correlated surfaces morphologies is studied using the Keldysh non-equilibrium Green function tech- 
nique. Due to the pseudo-periodicity of these structures, the energy-resolved transmission possesses 
pseudo-bands and pseudo-gaps. Channels with correlated surfaces exhibit wider pseudo-bands than 
their anti-correlated counterparts. By surveying channels with various combinations of material 
parameters, we found that a smaller transport mass increases the channel transmittivity and en- 
ergy bandwidth of the pseudo-bands. A larger quantization mass yields a larger transmittivity in 
channels with anti-correlated surfaces. For channels with correlated surfaces, the dependence of 
transmittivity on quantization mass is complicated by odd-to-even mode transitions. An enhanced 
threshold energy in the energy-resolved transmission can also be observed in the presence of sur- 
face roughness. The computed enhanced threshold energy was able to achieve agreement with the 
experimental data for Si(110) and Si(100) devices. 



I. INTRODUCTION 

In the literature, theoretical studies of the physics 
of surface roughness on electronic transport properties 
mainly focus on the linear response near thermodynamic 
equilibrium. In this regime, transport is diffusive and 
the electron dynamics is well described by a Boltzmann 
equation [T] or Kubo formula [2]. Once the pertubation 
Hamiltonian for surface roughness (Hsr) is formulated, 
the transition amplitude between electronic states can be 
computed through Fermi's golden rule. Surface rough- 
ness limited mobility can then be systematically calcu- 
lated. The theory on the form of Hsr traces back to the 
work by Prange and Nee [3] on magnetic surface states 
in metals. More recently, a systematic derivation of Hsr 
was discussed by Ando [4] in the context of electronic 
transport in semiconductors. 

It is well understood that this perturbation Hamilto- 
nian consists of two parts [5] : (i) local energy level fluctu- 
ations and (ii) local charge density fluctuations. The first 
term (i.e. local energy level fluctuations) is usually in- 
troduced phenomenologically [6] and explains the exper- 
imental observation that the surface roughness limited 
electron mobility in a quantum well scales with the film 
thickness (X5) and with the material quantization mass 
(m z ) according to T 6 -6 and ml respectively [6j [7J EJ E] • 
The latter term (i.e. local charge density fluctuations) is 
believed to be an important contribution to the degrada- 
tion of electron mobility in the high inversion charge den- 
sity regime [ID] . Although the treatment of the surface 
roughness problem is usually conducted in the framework 
of effective mass theory, a microscopic and self-consistent 
determination of Hsr can be obtained through density 
functional theory [11 . Another manifestation of surface 
roughness in quantum wells is the enhanced threshold 
energy, which has recently been observed experimentally 
[7] in quantum wells with thickness <4nm. These ex- 
periments show that the observed threshold energy does 



not follows the expected inverse quadratic scaling rela- 
tionship with T5. An objective of this paper is to explain 
why this deviation from quadratic scaling occurs. 

The physical effects of surface roughness on phase co- 
herent transport become very convoluted when the quan- 
tum well surfaces are roughened with random inhomoge- 
nuity of different scales [12 . We limit our study to 
phase coherent electronic transport through quantum 
wells with two special kinds of surface roughness mor- 
phology: (i) perfectly correlated surface roughness and 
(ii) perfectly anti-correlated surface roughness morpholo- 
gies. Phase coherent transport could be possible as de- 
vices are scaled into the deca-nanometer regimes [13] . In 
practise, one would expect a quantum well grown using 
the atomic layer deposition technique to produce a high 
degree of correlation between the two surfaces. Stud- 
ies of surface roughness scattering in the diffusive regime 
usually ignore such surface correlation effects [5]. Herein, 
we show that phase coherent transport through quantum 
wells with perfectly correlated or perfectly anti-correlated 
surfaces gives rise to distinctive features in the energy 
resolved transmission profile. The theoretical method 
we had employed is the Keldysh non-equilirium Green 
function (NEGF) approach [11 [15] [H Hi] within a fi- 
nite element and boundary element discretization scheme 

nanaED]. 

This paper is organised as follows. Section II discusses 
the NEGF formalism and methodology in a finite element 
discretization scheme. Section III examines the energy- 
resolved transmission characteristics for quantum films 
with correlated and anti-correlated surfaces. We discuss 
these results in comparison to the Kronig Penny model 
[2T] . Section IV studies the impact of quantization and 
transport masses on the transmission characteristics. Fi- 
nally, we compare our results to experimental values of 
the enhanced threshold energy for Si(100) and Si(110) 
devices. 



2 



II. THEORY AND MODEL 

The Landauer approach [22] [23] pictures a device 
within which dissipative processes is absent but coupled 
to perfect thermodynamic systems known as 'contacts'. 
This approach has been very successful in modeling phys- 
ical effects in myriad of problems in the field of meso- 
scopic physics [24]. When irreversible or energy dissi- 
pating processes are present in the device, a more so- 
phisticated quantum transport model such as the non- 
equilibrium Green function (NEGF) is needed to account 
for the coupling and transitions between the different 
quantum states in the system. The NEGF method was 
first formulated by Kadanoff and Baym [25] and Keldysh 
[T5] . Discussions of NEGF and its applications to con- 
densed matter phenomena can be found in textbooks by 
Datta [14], Haug [16] and Mahan [17] . 

This section summarizes the NEGF formalism applied 
to electronic transport through a 2D channel imple- 
mented using the finite element (FEM) and boundary 
element (BEM) discretization scheme [18] . Our choice 
of FEM over finite difference method is mainly because 
it can resolve the device's roughened surface geometry 
more efficiently with its flexible mesh. Our methods are 
similar to the ones developed by Havu et. al [20, 26] and 
Polizzi et al [19]. 
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FIG. 1: Illustration of the mesh used for simulation of lOnra 2D 
channel with (a) anti-correlated and (b) perfectly correlated surfaces. 
In our work, the roughness is characterized by only two parameters: 
amplitude (Ao) and wavelength (A) as depicted in (a). The meshes 
are generated with average distance of 0.25nra. 

Fig. [I] illustrates 2D channels with correlated and 
anti-correlated surface morphologies. The problem 
domain is denoted by with points represented by a 
2D spatial coordinate r=(x : z) G O. Vt is then partition 
into the interior domain ^ and exterior domains Q e j, 
where j=L, i?, 0. Q e L and Q e R denotes the left and 
right leads respectively, while fl e o denotes the remaining 
space. Each exterior domain Q e j shares the boundary 
with Qi denoted as dflij. The boundary of is simply 
dfti = YijdQij. The goal is to seek the numerical 
solution of the Green function in f^, denoted by G(r, r'). 
In our problem, the exterior domains Q e j consists of 
semi-infinite leads with known Green functions [T4] . 
Therefore, BEM can be applied to each of these exterior 



domain to account for its effect on the respective 
boundaries f^j . The purpose of FEM is to formulate the 
differential equation within £V 

Within the effective mass approximation, the 2D 
Hamiltonian that we are solving can be written as, 
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V r • [M(r)V r ^(r)l + V(r)^(r) 



(1) 



where r=(x,z) G and M(r) is a 2x2 effective mass 
tensor. In this work, we assumed an effective mass 



tensor with the following matrix elements; [M]n=m x 1 , 
[M]i2=0, [M]2i=0 and [M]22=^^ 1 - m x an d rn z are 
the transport and quantization masses respectively. The 
Green function of H is defined as, 



+ irj — G(r, r') — 5(r — r) : r £ 



n 



(2) 



where the boundary condition of outgoing waves is in- 
corporated by the introduction of r] — > + , i.e. G(r, r') 
is the retarded Green function. S(r — r') is the Dirac 
delta function. In the FEM scheme, we have the node- 
wise shape functions a^(r) as our basis functions [18] i.e. 
linear basis functions are employed in this work. Using 
the properties of a Dirac delta function of Eq. [2j we can 
write ah as follows, 

a*(r') = f ahWKW + e + ^G^r')^ 



:r'€f!, 



ai,(r)-V r • [M(r)V r G(r,r')] dQ 



(3) 

The second integral term on the R.H.S. of Eq. [3] contains 
a second order differential integrand. It can be reduced 
to first order via the identity, 

fiV • [/2V/3] = V • [/1/2V/3] - /2V/1 • V/ 3 (4) 

Eq. [3] then becomes, 

CLhir') = f a/,(r)(-y(r) + e + iry)G(r,r / )^ 

+ I (a h (r)^-M(r)V r G(r,r')) ■ nddSU 



/ 

J re 



JreSli 



M(r)V r a h (r) — • V r G(r,r)dQ 



(5) 



where ft in Eq. [5] is the normal vector to the boundary 
dfti. We note that in order to obtain V r G(r, r') along 
the boundary c^, one requires information about the 
Green function outside and within Qi. Recall that the 
Green function of the exterior domain Q e k, G e /c, has a 
simple analytical form [14]. For our purpose, we only 
need to know the explicit form of G e k along the boundary 
r ek edQ ik ; 



(~*ek(%eki ^eki %eki %ek) 
00 

= - Xm(Zek)Xm(z'ek) 



(6) 



m=l 



2sin(uj m x e k)exp(iuj m x e k) 
hv m , 
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where v m is the carrier velocity defined as Vm^hcOm/rrix 

and cu m =(l / h) T S /2m x (e — ft m ). Xm is the eigenstates of 
the confined modes in the leads corresponding to the 
eigenstate with energy K m . In order to incorporate the 
information of the exterior Green function into Eq. [5j we 
express the term V r G(r, r') in the last integral expression 
as [20], 



V r G(r, r) 



-L 



G{r e , 



„K 2 



M(r ej - ) V r V r . Gej (r ej , r)ddQ e 



: r G 



(7) 



With this replacement to Eq. [5j the calculus part of the 
problem is complete and we are ready to formulate the 



problem in matrix form. 

If we assume that the Green function can be expressed 
in terms of the FEM basis as follows; 



G(r,r') « ^(rKCOGi; 



(8) 



Then Eq. [5] can be formulated into a compact matrix 
equation by multiply both sides by a g (r f ) and integrate 
over r' . The matrix equation is, 



S = [eS - H - E] GS 



(9) 



where [S] g h= j a g (r f )ah(r f )dr f is commonly known as the 
overlap matrix. The explicit form for computing the ma- 



trix elements of E and H are given in Eq. 10 



[H] hi 



[/ 

[/ 



h 2 - 

<*h(r) (-V(r) + irj) a t (r) - —M(r)V r a h (r) • \/ r m(r)dQ 



h 2 - 

a h (r)—M(r) 



L 



h 2 

oti{r ek ) — M{r ek )VrVr ek G e k{r e k,r)ddVLek ) • nddQi 



(10) 



With all these quantities known, we are ready to com- 
pute the device observables through the Green function; 

G(e) = [eS-H-Z]- 1 (11) 

The device's local density-of-states is computed through 
the spectral function defined as, 

A(e) = C(e) 1 " [T L (e) + T R (e)] G(e) (12) 

where r=i[E — T,^] are the broadening functions to be 
computed individually for each leads [T4]. The diagonal 
elements of A(e) yields the local density-of-states. Fi- 
nally, the transmission is computed by taking the trace 
of the transmission function <l> defined to be [14], 

$(e) =r L (e)G(e)r R (e)G(ey (13) 

To facilitate our subsequent analysis, we unclustered the 
total transmission $ and examine only the transmission 
characteristics between mode m = 1,2 of left lead to 
mode n = 1,2 of right lead. This mode-to-mode trans- 
mission (3> mn ) is easily accomplished in our numerical 
scheme by noting that the lead self-energy can be un- 
folded into their respective modes: 

G e k = G\ k + G 2 ek + G\ k . . . (14) 

where G™ k is the self energy for mode m. This allows 
one to define a broadening function associated with each 
mode r™^™, from which the respective mode-to-mode 
transmission ^> mn can be computed. With this under- 
standing, we shall begin our numerical analysis. 

Figjl] illustrates a typical FEM mesh with (a) anti- 
correlated and (b) perfectly correlated surfaces used in 



our calculations. The FEM mesh is generated using 
the algorithm developed by Persson and Strang [27] 
based on the well-known Delaunay triangulation rou- 
tine. In our work, the roughness is characterized by si- 
nusoidal profiles with only two parameters: amplitude 
(Ao) and wavelength (A). These parameters are analo- 
gous to the root-mean-square roughness and the rough- 
ness auto-correlation length commonly employed in the 
literatures to describe the surface morphology [28]. In 
Fig. [2j we compare the energy-resolved transmission 
calculated with our FEM-BEM method to the scatter- 
ing matrix method (SMM) [29] using various mesh sizes 
d. In the SMM approach, the quantized energies are re- 
solved analytically unlike the FEM-BEM or finite differ- 
ence approaches. The simulations used a channel length 
I/G=5nm, an average film thickness Te=2nm, a trans- 
port mass m x =0.2 mo, and quantization mass m z =0.9mo. 
Anti-correlated surfaces described by roughness ampli- 
tudes t4o=0.1, 0.3, 0.5nm and wavelength A=2.5nm were 
considered. As shown in Fig. [2] the FEM-BEM re- 
sults are in satisfactory agreement with SMM. A very 
fine spatial discretization in the transport direction was 
employed for the SMM calculations. 



III. GENERAL FEATURES IN ENERGY 
RESOLVED TRANSMISSION 

The objective of this section is to perform a systematic 
analysis of the transmission characteristics of electronic 
transport through a 2D quantum channel with corre- 
lated and anti-correlated surfaces. As a model system, 
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FIG. 2: FEM and Scattering matrix method (SMM) simulation of 
a Lq = 5nm and Tb = 2nm channel with anti-correlated surfaces 
having roughness amplitudes Aq = 0.1, 0.3, 0.5nra and wavelength 
A = 2.5nra. We assumed a transport mass m x = 0.2rao and quanti- 
zation mass m z = 0.9rao i.e. based off Si material. SMM is a 'mode 
space' approach which the quantized energy is resolved analytically. 
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FIG. 3: Energy- resolved transmission (for&n) from FEM simulation 
of a Lq = lOnm and Tb = 2nm channel with various anti-correlated 
surfaces with roughness amplitudes Aq = 0.5nm and wavelength (a) 
A = 2.5nra, (b) A = 3.3nra and (c) A = 5.0nra. For each structure, 
we simulated a transport mass ofm x = 0.2rao (dash-dotted lines) and 
m x = O.lrao (solid lines) each with quantization mass m z = 0.9rao. 



we adopted a set of parameter values typical of 'end-of 
-road-map' devices [30] i.e. a Lc=10nm 2D channel with 
an average quantum film thickness of Te=2nm. For each 
class of devices, the effects of material properties and the 
degree of roughness on the energy resolved transmission 
characteristics were examined. The material bandstruc- 
ture was parameterized by a set of effective masses i.e. 
transport mass m x and quantization mass m z . Although 
microscopic details of the atomic structure of interfaces 
[TT] and related elemental defects (i.e. Si-Si and Si-O-Si 



bonds) [21] could be a source of interface scattering, it is 
not the main focus of this work. Ignoring these interfa- 
cial elemental defects, we study the geometrical effect of 
surface roughness on the transport characteristics. 



A. Transmission Through Anti-Correlated Surfaces 

Figj3] shows the energy-resolved transmission charac- 
teristics 3> mn =ii (first mode to first mode transmission) 
for devices with anti-correlated surfaces. Three differ- 
ent sets of roughness parameters were simulated: (a) 
Ao=0.5nm and A=2.5nm, (b) A =0.5nra and A=3.3nm, 
(c) Ao=0.bnm and A=5.0nm. For each set of rough- 
ness parameters, the following sets of material parame- 
ters were simulated: (dash-dotted lines) m z =0.9mo with 
m x =0.2m and (solid lines) m 2 =0.9m with m x =0.1m - 

Several general observations can be made about the en- 
ergy resolved transmission spectra: (i) there are regions 
of pseudo-gaps and pseudo-bands where transmissivity is 
relatively opaque and transparent respectively, (ii) within 
the pseudo-bands, there are camel-back structures which 
increases in number with increasing roughness frequency 
(i.e. decreasing A) and (hi) a delayed 'turn-on' of trans- 
mission called 'enhanced threshold', which increases with 
increasing roughness frequency The enhanced thresh- 
old energy is a geometrically derived property due to 2D 
quantization effects of the roughened morphology The 
enhanced threshold is zero for an unroughened quantum 
well channel. Note that the energy scale are referenced 
from the subband energy of the first mode in the lead, 
previously defined as K\. In addition, due to the symme- 
try of the problem, <£ mn 7^0 if and only if both m and n 
are odd/even numbers. 

Due to the anti-correlated surface morphology, the 
thickness of the 2D quantum film fluctuates from the 
source to drain contacts. Henceforth, we can visualise 
the electron as moving across the channel through an 
undulating energy landscape caused by the variable film 
thickness. This accounts for the appearance of camel- 
back structures in the pseudo-band region of the energy- 
resolved transmission as depicted in Fig (3] i.e. a signature 
of resonant tunneling. The undulating energy landscape 
(cqw(x)) can be modeled by a series of quantum wells 
(Eqw) as follows; 



£Qw(x) 



2m z [T B + 2A cos(2nx / 'A)]' 



= ^2e qw (x + jX) 

j 

(15) 



Eqw can be expanded as a Taylor series to give the fol- 
lowing leading order terms: 



Eqw(x) 



h 2 7T 2 



4h 2 n 4 A x 2 



2rn z (T B + 2A Y m z (T B + 2A ) 3 X 2 



U(x) 



(16) 



where U(x) is a unit pulse at —\/2<x<\/2. The sec- 
ond term in Eq. [TBlwith the kinetic operator will yield 
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FIG. 4: Available energy band states (indicated by shaded regions) for 
a Kronig Penny ID crystal with periodically varying square potential 
wells plotted with respect to the quantum well width A. The quantum 
well energy barrier is set to be OAeV. Effective mass of m x = 
0.15rao (solid) and m x = 0.2rao (patterned) are considered. 



the quantized energy levels of a harmonic oscillator with 
energies £ n =Huj(n + 0.5) where, 
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FIG. 5: Energy-resolved transmission (for <J>n and &12) from FEM 
simulation of a Lq = lOnm and Tb = 2nm channel with various per- 
fectly correlated surfaces with roughness amplitudes Aq = O.bnm and 
wavelength (a) A = 2.5nra, (b) A = 3.3nra and (c) A = 5. Own. For 
each structure, we simulated a transport mass of m x = 0.2rao (dash- 
dotted lines) and m x = O.lrao (solid lines) each with quantization 
mass m z = 0.9rao. 



8h 2 7i 4 A 



X 2 m x m z (T B +2A ) 



(17) 



From the above expression, we have c^ocra" - 5 . Indeed, 
FigJ3] reveals that a smaller m x will yield a wider en- 
ergy seperation between the peaks of the camel-back 
structure. This also translates to a larger pseudo-band 
bandwidth. For the structure illustrated in Figj3j the 
estimated first quantized energy using Eqjl7| is £i« 
0.1357eV, 0.1019eV and 0.0679eV for A=2.5nra, 3.33nra 
and 5nm, respectively. These estimates are in good 
agreement with threshold energies computed in Figj3j 
which is the energy needed for the first appearance of 
a transmission resonance. 

To explain the global features of pseudo-bands and 
pseudo-gaps in the energy-resolved transmission results, 
we shall first review the pertinent results from the 
Kronig-Penny model [51] . Consider a periodically vary- 
ing rectangular energy barrier as depicted in the inset 
of Fig. H] with barrier energy of Ub and a spatial period 
of A. The bandstructure (energy £ versus momentum n 
dispersion relation) of this ID periodic potential can be 
described by the following transcendental equation; 



±\tan 1 

A \ a 



(18) 



where a is given by; 



2cosh( — )cos( — ) 



(19) 



when 0<£<£4 and 



2cos( — )cos{ — ) 
k\ + k\ . ,k 2 \^ . ,ki\ s 



k\k 2 



-sin(—)sin( — ) 



(20) 



when ^>/7b. In addition, we have hki = ^/2m x \^\ and 
hk 2 = y / 2m x \£ - U b \. 

Fig. [4] surveys the bandstructure of a Kronig Penny 
ID crystal under different quantum well period A where 
regions with propagating states are shaded. The analy- 
sis set is conducted for transport masses of m x =0. 15mo 
(colored region) and m x =0.2mo (shaded region). The 
following observations can be made: (i) the energy band- 
width and threshold energy increase monotically with de- 
creasing quantum well period A, (ii) energy bandwidth 
decreases monotonically with decreasing quantum well 
period A and (iii) a smaller transport mass m x yields 
a larger energy bandwidth. Examination of the energy- 
resolved transmission in Figj3] shows that the pseudo- 
band and enhanced threshold in the quasi-periodic 2D 
structure with anti-correlated surface roughness exhibits 
a strikingly similar trend with the Kronig Penny model 
analysis. Based on these arguments, we conclude that 
the generic features of pseudo-band and pseudo-gap ob- 
served in the energy-resolved transmission of the 2D film 
with anti-correlated surfaces is a result of the film's quasi- 
periodicity. 
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FIG. 6: Intensity plot of local density of states due to the first mode 
of left contact i.e. G^T^G, each at different injection energies (mea- 
sured with respect to the lead's first quantized mode) as indicated. 
Device structure used is a Lq = lOnm and Tb = 2nm channel with 
anti-correlated surfaces with roughness amplitudes Aq = 0.5nm and 
wavelength X = 2.5nm. A transport mass ofm x = 0.2rao and quanti- 
zation massm z = 0.9rao are used. The corresponding energy-resolved 
transmission characteristics is plotted on the left. 
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the different symmetry in this case, <£ mn ^ for odd 
to even mode transition, unlike the situation for anti- 
correlated surfaces. Therefore, we plotted the transmis- 
sion characteristics 3>n and $12 in Fig. [5| although we 
will mainly focus on $n in the subsequent disussion. 

As shown in Fig. |5j the generic features of pseudo- 
band, pseudo-gap and enhanced threshold, induced by 
the quasi-periodicity are also present in structures with 
perfectly correlated surfaces. For each particular set of 
surface roughness parameters, the channel with perfectly 
correlated surfaces exhibits a distinctively larger pseudo- 
band than its anti-correlated counterpart. Since there is 
no variations of quantum well thickness along the chan- 
nel, the scattering in this case is purely a result of 2D 
geometrical and side- wall boundary effects. As a result, 
one would expect the carrier to feel a less undulating 
energy landscape in the perfectly correlated case. Ef- 
fectively, this translates to a smaller in the Kronig- 
Penny picture, which will then correpondingly yield a 
larger pseudo-band bandwidth. 

Fig. [6] and [7] are intensity plots for the local density- 
of-states {e)Y 1 L {e)G{e) due to injection of carriers from 
the source contact for the anti-correlated and corre- 
lated cases respectively. Brighter regions indicate higher 
density-of-states. The carriers are injected from the first 
mode eigenstates of the source lead. The surface rough- 
ness morphologies are both characterized by Ao=0.5nm 
and A=2.5nm. The energy resolved transmission spec- 
tra are plotted on their left. The local density-of-states 
at each of the resonance energies reveals localised high 
intensity patterns, which has its origin in interference ef- 
fects due to multiple scattering of waves. One observes 
that the number of localised spots increases with the in- 
dex of the resonance level. In the energy-resolved trans- 
mission, twice the number of resonance energy levels in 
the pseudo-band is observed in the perfectly correlated 
case as compared to anti-correlated case. This is at- 
tributed to the reduced symmetry of the channel with 
a perfectly-correlated surface morphology, which lacks 
the x-axis mirror symmetry of the anti-correlated sur- 
face. The degradation of the energy-resolved transmis- 
sion at higher energy is due to the appearance of the 
second mode at 0.3eV\ which opens up a transmission 
through $12. Therefore the resonance peak has maxi- 
mum transmission less than 1. 



FIG. 7: Similar to Fig^ except for a channel with perfectly-correlated 
surfaces. 



B. Transmission Through Perfectly- Correlated 
Surfaces 

Fig. [5] shows the energy-resolved transmission charac- 
teristic of a 2D channel with perfectly correlated surface 
roughness. The same sets of devices as performed in the 
anti-correlated case, with various surface roughness pa- 
rameters and transport mass m x , are simulated. Due to 



IV. EFFECTIVE MASSES ON 
TRANSMITTIVITY AND THRESHOLD ENERGY 



In this section, we study the impact of transport and 
quantization masses and surface roughness morpholo- 
gies on the general characteristics of the energy-resolved 
transmission through a deca-nanometer channel. Our 
analysis will be confined to the study of the first mode 
to first mode transmission, We begin by proposing 
a reasonable metric for the measure of the transmittivity 
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Roughness A mplitude(nm) Roughness A mplitude(nm) Energy wrt 1st Subband (eV) 



FIG. 8: Surverying of pseudo-band's bandwidth of the energy-resolved transmission for <E>n (see (a)) and the transmittivity (see (h)) for 
various transport masses (m x = O.lrao and 0.2rao, denoted by red and black lines respectively) and quantization masses (m z = 0.3rao, 0.5mo 
and 0.9rao, denoted by dash-dotted, dashed and solid lines respectively), plotted as a function of roughness amplitudes Aq. Simulated for a 
device with Lq = lOnm and an average Tb = 2nm channel with anti-correlated surface roughness morphology (X = 3.3nm). (c) depicts the 
energy-resolved transmission for similar devices but with Lq = 20nm (black) and 40mm (red), with various Ao=0.35nra and 0.50nra using 
effective masses of m x = 0.2rao and m z = 0.9rao. 
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FIG. 9: Effects of quantization mass (m z ) on the pseudo-band's 
bandwidth of the energy-resolved transmission for <E>n (see (a)) and 
the transmittivity (see (b) ) through a channel with perfectly-correlated 
surface roughness morphology plotted as a function of Aq. The device 
has Lq = 10nm and an average Tb = 2nm channel. Transport mass 
of m x = 0.2rao is used for all curves. The red curves in (b) denotes 
that higher modes transitions, i.e. <E>i2 and <E>22, are present within 
the pseudo-bands of these structures. 



K of a channel: 



K 



(21) 



where the pseudo-band bandwidth (BW) is defined as 
the energy difference between the last and first resonance 
peaks in the pseudo-band. Enhanced threshold energy 
is defined as the energy for the appearance of the first 
resonance peak with respect to the lead's first subband 
energy, i.e. h 2 7T 2 / (2m z Tg). This phenomenon has been 
observed experimentally [7J [9] and the suppression of its 
effect is pertinent to electronic device applications. Fig. 
|8jc) provides an illustration of the concept of BW and 
enhanced threshold energy. 

Fig. [8ja)(b) and |9ja)(b) surveys the BW and K 
for a deca-nanometer channel of (Te)=2nm with anti- 
correlated and perfectly correlated surface roughness 
morphologies respectively. Different sets of material pa- 
rameters, i.e. m X:Zl are employed in the study. The key 
results can be summarized as follows: (i) a smaller m x 
improves the pseudo-band's bandwidth (BW) and effec- 
tive transmittivity (K) for channels with anti-correlated 
and perfectly correlated surface roughness morphology, 
(ii) a larger m z improves the transmittivity for chan- 
nels with anti-correlated surface roughness but slightly 
degrades the transmittivity for channels with perfectly 
correlated surface roughness. The observation (i) is due 
to the increase of resonance linewidth as derived from the 
higher tunneling probability due to a smaller m x . 

In the diffusive regime, the surface roughness limited 
mobility for a quantum well scales proportionally with 
w m^A -2 , where A is the root-mean-square averaged 
fluctuation of the quantum well thickness [8 . In the 
channel with anti-correlated surfaces, A=\^Aq. As de- 
picted in Fig. [8jb), K increases with increasing m z . For 
the channel with perfectly correlated surface morphology, 
A is zero. Thus, the main source of scattering mechanism 
in the 'classical' sense is atrributed to a local fluctuation 
of wavefunction [5]. In general, a larger m z 'propagates' 
the electron closer to the surfaces and renders it more 
sensitive to the surface roughness morphology. This ex- 
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plains the small degradation of K with increasing m z . 
However, we must emphasize that the larger degrada- 
tion of K (red curves in^b)) for m z — 0.7mo,0.9mo is 
due to the appearance of the second mode, leading to a 
degradation of $11 while 3>i2 begins to increase. 

With a larger transmission bandwidth and relatively 
weak dependence of transmittivity on m zi channel with 
highly correlated surfaces is more optimal for electronic 
transport in the phase coherent regime. Especially for a 
channel with small quantization mass, i.e. III-V semicon- 
ductor alloys, the latter property is highly desireable for 
ballistic transport. Although these studies are conducted 
for a deca-nanometer channel, we expect the results to 
be consistent for longer channel. To confirm this proposi- 
tion, we considered the energy-resolved transmission for 
devices with Lc=20nm and 40nm in Fig. [8^c). The in- 
crease in channel length results in more resonance peaks 
(i.e. the number of peaks within the pseudo-band is pro- 
portional to the number of sinusoidal cycles in the rough- 
ness morphology of the channel) but with the global fea- 
tures of enhanced threshold and the pseudo-band's band- 
width intact. In particular, its energy-resolved transmis- 
sion exhibits similar enveloping characteristics for differ- 
ent Lc, where the envelope is characteristics of a given 
A . 

It has been reported that surface roughness will induce 
an additional threshold energy in quantum well devices 
and this effect had been systematically measured in ex- 
periments [9j [32] , i.e. enhanced threshold energy. There- 
fore, the observed threshold energy will be greater than 
that described by the more commonly understood body 
quantization effect according to h 2 7r 2 /(2m z T B ). Fig. 
To| compares the theoretically calculated threshold en- 
ergy of a roughened channel (i.e. perfectly correlated 
and anti-correlated surfaces morphology) with the avail- 
able experimental data for Si(llO) and Si(100) quantum 
well devices [7[ [9] . Assuming a roughness amplitude of 
Ao = 0.5nm and Ao = 0.6nm for the Si(llO) and Si(100) 
devices respectively, we are able to match the experimen- 
tal data for the range of quantum well thicknesses T#. A 
roughness wavelength of A = 2.5nm was assumed for all 
devices in our calculations. The corresponding threshold 
energy due to an unroughened channel are also shown 
(plotted as solid lines). A quantization mass (m z ) of 
0.9mo is employed for the Si(100) devices. For Si(llO), 
the effective mass tensor has off-diagonal terms in the 
direction normal to the quantum film surface. A uni- 
tary transformation is employed to decoupled them, as 
described by Stern and Howard [33] . Eventually, a quan- 
tization mass of 0.33mo is obtained for Si(llO) devices. 

The enhanced threshold energy is a geometrically de- 
rived property of 2D quantization effects. The 2D geom- 
etry of the roughened channel introduces additional lat- 
eral confinement which serves to enhance overall thresh- 
old energy. Fig. [6] and [7] illustrates this lateral confine- 
ment effect in channels with perfectly and anti-correlated 
surfaces. Note that these channels are constructed such 
that they retain the same volume as the unroughened 
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FIG. 10: Theoretically calculated threshold energy as a function of 
the averaged 2D quantum film's thickness compared with the exper- 
imental data for a 5/(110} and 5/(100) JTj devices. For our cal- 
culations, a roughness amplitude of Ao = 0.5nm and Ao = 0.6nm 
are able to describe the experimental data for 5/(110) and 5/(100) 
devices respectively where a roughness wavelength of A = 2.5nm 
was assumed for both cases. The threshold energy for a channel with 
unroughened (solid lines), anti-correlated (dashed lines) and perfectly 
correlated (dotted lines) surfaces are plotted. 



channel. Eq. [TT] describes the approximate enhanced 
threshold energy due to surface roughness for the case of 
channels with anti-correlated surface morphology. Eq. [17 
tells us that the enhanced threshold energy is related to 
the material parameters according to (m x m 2 ) -1 / 2 . Fur- 
thermore, contrary to the T^ 2 dependency in the usual 
case of body quantization, the enhanced threshold en- 
ergy exhibits a T B 3 ^ 2 dependency as shown in Eq. [l7| 



From the viewpoint of device performances, this trans- 
lates to an additional threshold voltage shift. A small 
device threshold voltage shift is an important criterion 
to suppress the on-chip device-to-device electrical prop- 
erties variations l34l. 



V. CONCLUSION 

The study of the impact of surface roughness on 
electronic transport had usually focused on the dissi- 
pative regime, with transport dynamics well governed 
by the classical Boltzmann's transport equation. Elec- 
tronic transport through a channel with roughened 
surfaces in the phase coherent regime is less under- 
stood. Using the Keldysh non-equilibrium Green func- 
tion approach within a FEM-BEM numerical scheme, 
we performed a systematic study of quantum transport 
through a deca-nanomater length quantum well chan- 
nel with perfectly correlated and anti-correlated sur- 
faces. Due to the pseudo-periodicity in these simulated 
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structures, their energy-resolved transmission possesses 
pseudo-bands, pseudo-gaps and an enhanced threshold 
energy. Channels with perfectly correlated surfaces ex- 
hibit wider pseudo-bands than their anti-correlated coun- 
terparts. Perfectly correlated channels also permit odd- 
to-even mode transition, which is not allowable in chan- 
nel with anti-correlated surfaces. 

An effective transmittivity in these structures is de- 
rived, by computing the average transmission over the 
range of energy within the pseudo-band. By surveying 
channels with various material parameters combinations 
(i.e. m x and m z ), we found that a smaller transport mass 
m x is beneficial for the transmittivity of the channel and 
serves to increase the energy bandwidth of the pseudo- 
band. The observation of the contrasting trends in the 
dependence of transmittivity on m z for anti-correlated 
and perfectly correlated surfaces is interesting. A quan- 
tum well with perfectly correlated surface is more opti- 
mal for channel material with smaller m z . Technically 
speaking, a sufficiently correlated surface can be engi- 
neered via techniques like atomic layer deposition. Recall 
that the 'classical' perturbing Hamiltonian due to surface 
roughness can be attributed to a local energy level fluc- 
tuation and a local fluctuation of charge density [5 . On 
a general note, one could then say that channels with 
anti-correlated surfaces emphasize the former scattering 



mechanism, while the perfectly correlated surfaces em- 
phasize the latter mechsmism. 

Lastly, we studied the phenomenon of enhanced 
threshold voltage shifts. Excellent corroboration with 
the experimental data was obtained. Enhanced thresh- 
old voltage shifts exhibits a T B 3 ^ 2 dependency and its 
contribution to the total threshold energy is significant 
in the small 7# regime. Therefore, suppression of the 
device-to-device threshold voltage variations in quantum 
well channels with small quantization mass will present 
considerable challenge for the semiconductor device in- 
dustry [30] . 
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